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ANALYSIS  OF  THE  ORBIT  OF  1966-63A  NEAR  1 5TH-ORDER  RESONANCE 

by 

Doreen  M.  C.  Walker 
A.  N.  Winterbottom 


SUMMARY 

The  orbit  of  1966-63A  has  been  determined  at  32  epochs  between 
October  1976  and  April  1977,  using  the  RAE  orbit  refinement  program  PROP 
and  some  1500  US  Navy  observations. 

The  orbit  passed  through  exact  15th-order  resonance  on  1976  December  8 
and  the  32  values  of  inclination  have  been  analysed  to  obtain  lumped  15th- 
order  geopotential  coefficients, 

<3-0,1  9_0,1 

10  C,  =  36900  ±  9700  and  10  S.,  =  12200  ±  1700  . 

lb  i  b 

When  these  values  are  used  (with  22  other  pairs  of  values)  to  solve  for 
individual  harmonic  coefficients,  there  is  a  definite  improvement  over  any 
previous  solution. 

It  was  not  possible  to  successfully  analyse  the  values  of  eccentricity 
as  the  effects  of  solar  radiation  pressure  were  too  great  and  the  actual 
change  due  to  resonance  very  small. 
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INTRODUCTION 


The  satellite  I966-63A,  otherwise  known  as  orbiting  vehicle  1-8  (OVI-8),  was  launched 
on  14  July  1966.  It  was  a  light  'spherical  skeleton',  mass  10.4  kg,  and  had  a  diameter  of 
9.14  m.  The  initial  apogee  and  perigee  heights  were  1013  km  and  998  km,  and  the  inclina¬ 
tion  was  144.27°.  The  satellite  was  in  orbit  for  llj  years,  decaying  on  4  January  1978  *. 

Jkhe  orbit  of  1 966-6  3A  had  an  initial  period  of  105.12  minutes  and  slowly  contracted 
under  the  influence  of  drag  to  pass  through  !5th-order  resonance  in  December  1976.  A 
satellite  experiences  15th-order  resonance  when  the  track  of  the  satellite  over  the  Earth 
repeats  aftei  i5  revolutions,  and  if  the  decay  rate  is  slow  enough,  an  orbit  passing 
through  the  resonance  can  be  appreciably  perturbed  by  the  effects  of  !5th-order  harmonics 
in  the  geopotential.  Recently  23  resonant  orbits,  at  various  inclinations,  have  been 

_ 

analysed  to  determine  harmonic  coefficients  of  order  15  and  degree  15,  16,  17,  ...  35 
The  satellite  1966-63A  was  used  in  this  analysis,  the  lumped  harmonic  coefficients  being 
obtained  from  the  analysis  of  25  values  of  inclination  from  US  Navy  orbits  kindly  supplied 
by  the  Naval  Research  Laboratory.  The  fit  of  the  theoretical  curve  to  the  values  of 
inclination  was  not  too  satisfactory,  sea_Eig  48  of  Ref  2,  and  the  analysis  of  the  eccen¬ 
tricity  was  not  attempted  as  the  values  were  not  considered  accurate  enough.  The  aim  of 
this  Report  is  to  compute  more  accurate  and  if  possible  more  frequent  orbits  over  the 

resonant  period  using  the  raw  observations  supplied  by  the  US  Navy  and  the  RAE  orbit 

■  -« 

refinement  program  PROP  in  the  PR0P6  version  .  — - 

The  orbits  determined  should  yield  more  accurate  lumped  coefficients  from  the 
analysis  of  inclination  and  possibly  allow  an  analysis  of  the  eccentricity.  This  satellite 

will  also  serve  to  test  the  PROP  program  in  determining  orbits  of  higher  inclination.  The 

»  o  A 

highest  so  far  determined  using  PROP  has  been  98.68  .  Finally  it  may  also  test  the" 

accuracy  of  the  method  for  restoring  the  odd  zonal  harmonic  oscillation  which  is  removed 

by  the  US  Navy  from  their  values  of  eccentricity3’^,  and  so  help  future  analysis. 

2  THE  OBSERVATIONS  AND  ORBITS 

The  orbit  of  1966-63A  has  been  determined  at  32  epochs  from  1481  observations 
between  3  October  1976  and  29  April  1977.  The  observations  used  were  Navspasur  observa¬ 
tions  supplied  by  the  Naval  Research  Laboratory,  and  they  have  a  topocentric  accuracy  of 
about  2  minutes  of  arc.  These  were  the  only  observations  available,  as  this  satellite 
was  not  placed  on  the  priority  list  for  visual  observers,  because  its  inclination,  144°, 
meant  that  it  only  reached  36°  latitude  and  was  therefore  not  easily  visible  from  the 
United  Kingdom. 

The  orbits  were  determined  at  approximately  weekly  intervals  and  the  orbital  elements 
at  each  of  the  32  epochs  are  listed  i.>  Table  I,  with  the  standard  deviations  below  each 
value.  The  epoch  for  each  orbit  is  at  00  hours  on  the  day  indicated,  and  the  PROP  program 
fits  the  mean  anomaly  M  by  a  polynomial  of  the  form 


M„  ♦  M.t  +  M-t2  +  M.t3  +  M.t4  +  Mct5 
0  12  i  4  3 


Table  X  (concluded) 


f 
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where  t  is  the  tine  measured  from  epoch  and  the  number  of  M-coefficients  used  depends  on 
the  drag.  For  1966-63A,  which  was  in  a  nearly  circular  orbit  at  a  height  of  about  600  km, 
Mq-M2  were  the  only  coefficients  required  on  14  orbits;  another  14  orbits  required 
to  be  added;  and  four  orbits  needed  M^-M^  .  The  value  of  e  ,  the  parameter  indicating 
the  measure  of  fit,  was  between  0.31  and  0.98,  having  an  average  value  of  0.50.  The  obser¬ 
vations  were  therefore  fitted  in  a  satisfactory  manner  on  all  the  orbits. 

The  values  of  inclination,  i  ,  had  standard  deviations  varying  between  0.0005°  and 
0.0013°,  the  average  being  0.0007°  corresponding  to  about  80  m  in  distance.  The  standard 
deviations  for  eccentricity,  e  ,  ranged  between  0.00001  and  0.00004,  the  average  value 
being  0.00002,  equivalent  to  130  m  in  distance. 


The  values  of  inclination  from  the  32  epochs  in  Table  1  have  been  plotted  as  circles 
in  Fig  I  and  a  connecting  line  has  been  drawn  through  the  points.  The  values  of  inclina¬ 
tion  from  the  US  Navy  orbits  over  the  same  period  are  plotted  as  crosses.  The  two  sets  of 
values  show  the  same  trends,  the  PROF  values  being  slightly  higher  than  the  US  Navy  values. 
Only  two  of  the  US  Navy  values  differ  from  the  line  joining  the  PROP  values  by  more  than 
the  sum  of  their  standard  deviations,  if  the  usual  value  of  0.003°  is  taken  for  the 
standard  deviation  of  the  US  Navy  values.  This  result  is  most  satisfactory  and  confirms 
the  use  of  0.003°  as  the  standard  deviation  for  the  US  Navy  values  of  inclination.  At  this 
inclination  the  US  Navy  values  appear  to  have  a  bias  of  0.0018°  relative  to  PROP. 


3  I 5TH-ORDER  RESONANCE 

3. I  Analysis  of  inclination 

The  theoretical  equation  for  the  variation  of  inclination  at  )5th-order  resonance 

,  2 
may  be  written  as 


di 

dt 


-M-)15  r< 

jm  i  Va / 


(15  -  cos  i)F 


J>,1 

15,15, 7/C15 


sin  *  -  S 


.0,1 


15 


cos  4 


I 

1 


+  terms  in 


(Ue)lql 

(111)’ 


(y*  -  qio) 


0- 


(2) 


where  t  the  resonance  angle  is  given  by 


♦  ■  u  ♦  M  ♦  15(0  -  v) 


(3) 


and  at  exact  resonance  $»0  .  The  variations  of  4  and  $  are  given  in  Fig  2,  and  the 
definition  of  all  the  symbols  used  in  equations  (2)  and  (3)  is  the  same  as  in  Ref  2.  In 

equation  (2)  only  the  main  term  is  given  explicitly:  the  terms  involving  sine  and  cosine 

2*  will  probably  not  be  required,  as  they  are  associated  with  harmonics  of  order  30  and 

should  be  considerably  smaller  than  those  of  order  15;  and  with  e  <  0.003  during  the 
resonance  period,  the  e-terms  will  be  small  also. 

Hie  32  values  of  inclination  in  Table  I  were  cleared  of  zonal  harmonic  and  luni- 
solar  perturbations  using  the  PROD^  computer  program  with  one-day  integration  steps,  and 
the  2  perturbations  were  also  removed.  The  resulting  values  of  inclination,  with 
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the  standard  deviations  quoted  in  Table  I,  were  then  fitted  with  equation  (2)  in  integrated 

g 

'.orm  by  the  THROE  computer  program,  which  removes  perturbations  due  to  atmospheric  rota¬ 
tion  and  lunisolar  precession  of  the  Earth's  axis.  The  density  scale  height  H  was  taken 

9 

as  78  km,  appropriate  to  a  height  of  624  km,  0.21H  above  perigee  ,  and  the  atmospheric 
rotation  rate,  A  was  taken  as  0.8  rev/day. 

This  first  run  by  THROE  using  just  the  (y,q)  =  (1,0)  terms  as  given  explicitly  in 
equation  (2)  produced  a  rather  poor  fit  to  the  values  of  inclination.  The  values  of 
lumped  harmonics  obtained  were: 

q_0,l  Q.0.1 

10  C[5  -  86700  ±  15100  10  S|5  =  -  7900  ±  7300  , 


with  e  «  2.1  .  The  values  of  inclination  cleared  of  the  perturbations  already  mentioned 

are  plotted  in  Fig  3  together  with  the  THROE  fitting  with  (y,q)  «  (1,0).  From  Ref  2  it 

9  0,1  9*0,1 

was  known  that  the  values  of  10  and  10  5|,_  should  be  of  order  49000  and  19000 

respectively.  Further  runs  were  tried  with  more  coefficients  and  with  selected  relaxations 

.0,1 

in  accuracy;  but  the  extra  coefficients  were  not  determined,  the  lumped  coefficients  C,, 
.0,1  15 
and  Sj,.  were  no  nearer  to  the  values  required  and  there  still  remained  a  variation  in 

inclination  which  had  not  been  successfully  removed  or  fitted. 


At  this  point  it  was  realized  that  as  1966-63A  was  a  spherical  skeleton,  it  would 
probably  be  affected  by  solar  radiation  pressure.  Dr  P.  Moore  at  the  Earth  Satellite 
Research  Unit  of  the  University  of  Aston  in  Birmingham  has  developed  a  program  for  comput¬ 
ing  the  effects  of  solar  radiation  pressure  on  the  orbital  elements,  and  kindly  ran  this 
computer  program  for  us,  so  that  the  values  of  inclination  could  then  be  cleared  of  this 
further  perturbation  and  the  THROE  fitting  tried  again. 


In  order  to  obtain  the  correct  change  in  inclination  due  to  solar  radiation  pressure 

a  value  of  the  mass/area  ratio  of  the  satellite  has  to  be  supplied.  Here  it  was  found 

2 

that  a  mass/area  ratio  of  2  kg/m  resulted  in  the  best  THROE  fitting  to  the  values  of 
inclination.  The  values  of  the  coefficients  were: 


q_°»'  q.0,1 

10’c)5  -  36900  ±  9700  10  S15  *  12200  1  4200  . 


with  e  "  1.3  .  This  final  fitting  and  the  values  of  inclination  with  all  known  perturba¬ 
tions  removed  are  plotted  in  Fig  4.  In  this  fitting  the  standard  deviations  of  the  values 

at  MJD  43134  and  43140  were  increased  by  factors  of  3  and  2  respectively. 

_0,1 

When  the  new  value  of  C.,  is  used  in  place  of  the  old  value  for  I966-63A  in  the 

5  2 

solution  for  individual  15th-order  harmonics  ,  there  is  an  improvement  of  about  IZ  in  the 

„0,l 

fit.  When  the  new  value  is  used,  however,  the  improvement  is  more  impressive:  the 

previous  value,  -6000  ±  7600,  had  to  have  its  standard  deviation  increased  by  a  factor  of 

.0,1 

4  m  order  to  keep  its  weighted  residual  below  1.5,  whereas  the  new  value  of  3|j  with 

its  standard  deviation  unchanged  gives  a  weighted  residual  of  1.39.  In  the  new  solutions 

for  individual  15th-order  harmonics,  the  values  do  not  differ  from  those  previously 
2 

obtained  by  more  than  1/10  of  the  standard  deviations.  (The  new  values  are  to  be 
regarded  as  better  than  the  old  ones:  but  they  will  not  be  quoted  here,  to  avoid  the 
confusion  that  would  be  caused  by  a  proliferation  of  nearly  identical  solutions.) 
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3.2  Analysis  of  eccentricity 

The  perturbations  in  the  eccentricity  due  to  zonal  harmonics  and  atmospheric  drag 
must  be  removed  before  an  attempt  is  made  to  analyse  the  resonance  effect.  The  air  drag 
model  within  the  THROE  program  takes  no  account  of  the  day-to-night  variation  in  density, 
so  a  correction  was  applied  to  the  values  of  eccentricity  using  the  method  detailed  in 
Ref  11.  As  solar  radiation  pressure  had  affected  the  values  of  inclination,  its  effect 
on  eccentricity  would  be  expected  to  be  even  more  important,  so  these  perturbations  were 

also  removed  using  the  data  supplied  by  Dr  P.  Moore  of  Aston  University,  with  the  mass / 

.  2 
area  ratio  again  assumed  to  be  2  kg/m  . 

The  values  of  eccentricity  cleared  of  all  known  perturbations  except  those  due  to 
resonance  are  plotted  as  circles  in  Fig  5.  These  values  of  eccentricity  were  then  fitted 
with  equation  (12)  of  Ref  2  in  integrated  form  using  THROE,  taking  (y,q)  “  (1,1)  and 
(1,-1)  terms.  Tl.  fitted  curve  given  by  this  THROE  run  is  shown  as  a  full  line  in  Fig  5. 
This  curve  did  not  fit  the  plotted  points  at  all  well  (e  =  11.4)  and  it  was  evident  that 
there  was  still  a  perturbation  in  the  values  of  eccentricity  unaccounted  for.  This  is 
probably  because  the  solar  radiation  pressure  perturbation  is  more  complex  than  has  been 
assumed  and  is  not  adequately  removed  by  taking  a  constant  mass/area  ratio. 

A  further  THROE  run,  non-fitting,  with  the  values  of  lumped  harmonic  coefficients 
calculated  from  the  values  of  individual  coefficients  of  Ref  2,  gave  the  variations  shown 
by  the  broken  curve  in  Fig  5.  This  curve  indicates  that  resonance  causes  very  little 
change  in  eccentricity,  corresponding  to  a  maximum  change  in  distance  of  around  400  m, 
and  therefore  has  a  much  smaller  effect  than  solar  radiation  pressure  (4  km),  zonal  harmo¬ 
nics  (7  km),  or  air  drag  (basic  4  km,  day-to-night  1.5  km).  In  view  of  the  uncertainty 
of  the  solar  radiation  pressure  perturbation,  there  seems  little  hope  of  isolating  the 
resonance  variation,  which  is  only  three  times  the  average  standard  deviation  of  the 
values.  So  the  attempt  to  analyse  the  eccentricity  was  abandoned. 

4  COMPARISON  OF  PROP  AND  US  NAVY  VALUES  OF  ECCENTRICITY 

The  values  of  eccentricity  in  Table  I  have  been  plotted  as  circles  in  Fig  6  after 
the  perturbations  due  to  zonal  harmonics  had  been  removed  using  the  PROD^  computer  program 
with  one-day  integration  steps.  These  values  of  eccentricity  are  joined  with  a  connecting 
curve.  With  the  US  Navy  values  of  eccentricity  we  first  have  to  restore  that  part  of  the 
odd  zonal  harmonic  oscillation  which  the  US  Navy  remove.  This  is  done  by  using  the 
ELTRAN  2  program,  and  should  render  the  US  Navy  values  equivalent  to  the  PROP  values. 

These  adjusted  US  Navy  values  of  eccentricity  were  then  cleared  of  perturbations  due  to 
zonal  harmonics  in  the  same  way  as  the  PROP  values  and  plotted  as  crosses  in  Fig  6. 

In  previous  analyses  of  US  Navy  values  of  eccentricity  there  has  often  been  a 
residual  variation  correlated  with  w  ,  particularly  for  satellites  at  inclinations 
greater  than  50°.  However,  Fig  6  shows  no  sign  of  such  an  oscillation,  which  would  have 
a  period  of  about  45  days.  Also,  the  general  agreement  between  the  US  Navy  and  PROP 
values  is  fairly  satisfactory:  if  the  standard  deviations  are  taken  as  0.00004  and 
0.00002  respectively,  16  of  the  25  US  Navy  values  differ  from  the  curve  through  the  PROP 
values  by  less  than  the  sum  of  the  standard  deviations.  The  inclination  of  I966-63A 
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(144.16°)  is  equivalent  to  35. 84°,  so  Fig  6  reinforces  the  conclusion  that  deficiencies 
in  the  conversions  made  by  ELTRAN  2  are  only  apparent  for  orbits  inclined  at  more  than  50° 
to  the  equator. 

Fig  6  also  serves  to  show  the  very  large  variation  in  eccentricity  even  after  the 
removal  of  zonal  harmonic  perturbations:  the  large  increase  at  the  end,  amounting  to  about 
9  km,  is  probably  due  to  solar  radiation  pressure,  and  this  is  why  we  failed  to  fit  the 
eccentricity  variation  in  section  3.2,  where  a  constant  mass/area  ratio  was  assumed,  which 
gave  only  a  4  km  change.  It  seems  that  the  mass/area  must  have  decreased  in  March  1977. 

5  CONCLUSIONS 

The  orbit  of  I966-63A  has  been  determined  at  32  epochs,  between  October  1976  and 
April  1977,  when  the  effects  of  resonance  were  being  felt.  About  1500  observations  made 
by  Navspasur  and  kindly  supplied  by  the  Naval  Research  Laboratory  were  used  in  the  deter¬ 
mination.  The  32  orbits  obtained  are  given  in  Table  1 ,  the  values  of  inclination  and 
eccentricity  having  average  standard  deviations  equivalent  to  80  m  and  130  m  in  distance 
respectively. 

This  was  the  first  time  the  PROP  program  had  been  used  to  determine  orbits  with 
inclination  greater  than  100°,  The  program  ran  smoothly  and  no  problems  were  experienced. 

The  32  values  of  inclination,  cleared  of  all  known  perturbations  except  those  due 
to  resonance,  were  fitted  with  a  theoretical  curve  and  the  following  values  of  lumped 
15th-order  coefficients  were  obtained: 

„_°,1  o.0,1 

10  C|5  -  36900  ±  9700  10  si5  “  12200  ±  4700  • 


As  I966-63A  was  a  spherical  skeleton,  it  was  found  necessary,  for  the  first  time  in  a 

resonance  analysis,  to  remove  perturbations  due  to  solar  radiation  pressure  before  the 

theoretical  curve  could  be  fitted.  This  perturbation  was  removed  by  utilising  a  program 

developed  by  Dr  P.  Moore  at  Aston  University:  the  subsequent  fit  gave  very  satisfactory 

-0,1  ,0,1 

results,  and  the  lumped  values,  and  S,,.  yielded  better  solutions  for  individual 


15 


15th-order  harmonics,  especially  through  the  new  value  of  S 


.0,1 

15 


The  analysis  of  the  values  of  eccentricity  was  disappointing,  because  solar  radia¬ 
tion  pressure  had  such  a  strong  effect  that  a  determination  of  the  lumped  harmonics  was 
impossible.  However,  the  PROP  values  of  eccentricity  were  useful  in  showing  that  the 
US  Navy  values  of  eccentricity,  after  conversion  with  the  ELTRAN  2  program,  were  consist¬ 
ent  with  the  PROP  values. 
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Fig  3  Values  of  inclination,  cleared  of  zonal  harmonic,  hinisolar  and  atmospheric 
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with  fitted  curve 
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Fig  8  Values  of  eccentricity  from  PROP  and  US  Navy  orbits,  with  zonal  harmonic 
perturbations  removed 
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